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GENERAL SOLUTIONS OF SUMS OF CONSECUTIVE CUBED 
INTEGERS EQUAL TO SQUARED INTEGERS 


VLADIMIR PLETSER. 


Abstract. All integer solutions (M, a , c) to the problem of the sums of M 
consecutive cubed integers (a + i) 3 (a> 1, 0 < i < M — 1) equaling squared 
integers c 2 are found by decomposing the product of the difference and sum 
of the triangular numbers of (a + M — 1) and (a — 1) in the product of their 
greatest common divisor g and remaining square factors S 2 and a 2 , yielding 
c = gSa. Further, the condition that g must be integer for several particular 
and general cases yield generalized Pell equations whose solutions allow to find 
all integer solutions (M, a, c) showing that these solutions appear recurrently. 
In particular, it is found that there always exist at least one solution for the 
cases of all odd values of M, of all odd integer square values of a, and of all 
even values of M equal to twice an integer square. 
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1. Introduction 


It is known since long that the sum of M consecutive cubed positive integers starting 
from 1 equals the square of the sum of the M consecutive integers, which itself equals 
the triangular number Am of the number of terms M, 

2 

'M(M + 1)' 


(i.i) 



_ A2 
— 


for VM £ Z + . The question whether this remarkable result can be extended to 
other integer values of the starting point, i.e. whether the sum of consecutive 
cubed positive integers starting from a ^ 1 is also a perfect square 


M -1 

( 1 . 2 ) J2(a + if = c 2 

i=0 

has been addressed by several authors but has received so far only partial answers. 
With the notation of this paper, Lucas stated [15] that the only solutions for M = 5 
are a = 0,1,96 and 118 (missing the solution a = 25, see further Table 1), and that 
there are no other solutions for M = 2 than a = 1. Aubry showed [1] that a 
solution for M = 3 is a = 23, c = 204, correcting Lucas’ statement that there 
are no other solution for M = 3 than a = 1. Other historical accounts can be 
found in [7]. Cassels proved [3] by using the method of finding all integral points 
on a given curve of genus 1 y 2 = 3x ( x 2 + 2) (with x = a + 1, y = c in this paper 
notations), that the only solutions for M = 3 are a = —1,0,1 and 23. Stroeker 
obtained [37] complete solutions for 2 < M < 50 and M = 98, using estimates of 


l 
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lower bound of linear forms in elliptic logarithms to solve elliptic curve equations 
of the form Y 2 = X 3 + dX where d = n 2 (n 2 — l) /4, X = nx + n (n — 1) /2, 
Y = ny (with n = M, x = a, y = c in this paper notations). The method 
reported, although powerful, appears long and difficult and caused some problems 
for the cases M = 41 and 44. Stroeker remarked also that M = a = 33 with 
c = 2079 = 33 x 63. This is not the single occurrence of M = a, as it occurs also 
for M = a = 2017,124993, 7747521,... (see [24, 25]). 

One of the reasons that these previous attempts to find all solutions (M, a, c) to 
(1.2) were only partially successful was most likely due to the approach taken to 
start the search for solutions for single values of M, one by one and in an increasing 
order of M values. The method proposed in this paper is to tackle the problem in 
a different way and instead of looking at each individual values of M one by one, 
to consider the problem in a more global approach by comparing different sets of 
known solutions, and instead of listing solutions in increasing order of M values, 
to look at two other parameters, 6 and a. defined further. This new beginning 
then leads to a more classical approach using general solutions of Pell equations, 
that allows to find all solutions in (M, a, c) of (1.2) for all possible cases. Note 
that Pell equations were already used previously by various authors (e.g. Catalan 
[4, 5], Cantor [2], Richaud [34]) in the 19th century in attempts to solve the present 
problem, albeit without reaching a complete resolution of the problem. 

The approach proposed in this paper includes three steps. 

Step 1 in Section 2 is based on the decomposition of the product of the difference A 
and the sum E of the two triangular numbers of (a + M — 1) and (a — 1) in simple 
factors, g, 5, a where g = gcd (A, E), allowing to find general expressions of a and c 
in function of M, S and er that are always solutions of (1.2). 

In step 2 in Section 3, some conditions on <)' and er are explored to obtain three 
particular cases of solutions yielding specific expressions of M in function of k, 
Vfc £ Z + , including the case of M taking all odd positive integer values. Section 4 
recalls some basics on Pell equation solutions to introduce the third step. 

In step 3 in Section 5, a general solution in (M, a, c ) is found for all values of 6 and 
er, based on solutions of simple and generalized Pell equations involving Chebyshev 
polynomials, allowing to find all the solutions in (M, a, c ) to (1.2). As an alternative 
to step 3, recurrence relations are deduced in Section 6. Section 7 summarizes all 
findings. 

Trivial solutions are not considered here. For instance, M < 0 is meaningless; for 
M = 1, the only solution is obviously a = a 2 , yielding (M, a,c) = (l,a 2 ,a 3 ); 
therefore, we consider only M > 1. If a < 0, there are no solutions if M < (1 — 2a); 
if M = (1 — 2a), the only solution is (M, a, c) = ((1 — 2a), a, 0); if M > 1 — 2a, 

the solutions are (M, a, c) = ^M, (1 — a), 2 M+a _i — A 2 _ a ^, i.e. a shift of the 
first term from a negative value a to a positive value (1 — a) with a reduction of 
the number of terms (M + 2a — 1). If a = 0 or 1, the classical solutions (1.1) 
are (M, a, c) = (M, 0, A| f _ 1 ) or (M, 1,A|^). Therefore, we limit our search to 
solutions for M > 1 and a > 1. 
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2. Step 1: A General theorem 


Theorem 1. For\/5 £ Z + , 3a, a, M, c £ Z + , k £ Q + , such as 5 < a, gcd (S, a) = 1, 
k = (cr/<5) > 1, 


( 2 . 1 ) 
holds if 
( 2 . 2 ) 
(2.3) 


M—l 

E (°+*) 3 = c 2 

i—0 


M (k 2 ~ 1) + 1 + JM 2 (k 4 - 1) + 1 
° ~ 2 

c = ^k 2 M + \JM 2 (k 4 — 1) + 1^ 


Proof. For a, M, c, A, S, g, A', S’, 6, a, C £ Z + , i,k £ Z*, k £ Q + , with M > 1, 
0 < i < M — 1, k = (cr/J) > 1, the sum of cubes of M consecutive integers (a + *) 
for i = 0 to M — 1 can be written successively as 


M -1 

( 2 - 4 ) + 

i=0 

(2.5) 

( 2 . 6 ) 

(2.7) 


a-\-M— 1 a —1 

E 

i —0 i=0 

A 2 — A 2 

(A a+ M-i — A a _i) (A a+ M-i + A a _i) 

(A) (£) 


where A = A a+ M-i — A a _i and £ = A a+ M-i + A a _i, i.e. the difference and the 
sum of the triangular numbers of (a + M — 1) and (a — 1), with obviously A < £, 
that can also be written 


( 2 . 8 ) 

(2.9) 


A 

£ 


M 



M - 1 \ 
2 ) 


a 2 + a (M - 1) + 


M (M - 1) 
2 


Let g = gcd (A, £), yielding A = g A' and £ = <?£'. For (2.1) to hold, c 2 = g 2 A'£' 
and as A' and £' are coprimes and their product must be square, both must be 
integer squares, i.e. A' = 6 2 and £' = cr 2 , with gcd (6, a) = 1 and S < a, yielding 


( 2 . 10 ) 


c = gSa 


From (2.8) and (2.9), one has then respectively 


( 2 . 11 ) 

( 2 . 12 ) 


2a + M - 1 


2 9 S 2 
M 


(M 2 - 1) 


where the + sign is taken in front of the square root in (2.12) as 2 a + M > 1. 
Solving for g yields then 

(2.13) g = H (k 2 M + \J M 2 (k 4 - 1) + l) 

Replacing in (2.11) or (2.12) yields then directly (2.2) and in (2.10) yields directly 
(2.3). ' ” □ 
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3. Step 2: Three Particular Solutions and a Generalization 

For which values of M does (2.1) hold ? Answers can be found for at least three 
particular cases. Other general approaches are developed further. 

Table 1 gives for 1 < M < 45 and 1 < a < 10 5 , the first values of M,a and 
associated A, E,g,<5, a and c values. For M and a < 10 5 , there are 892 solutions 
(M, a, c), given in [26, ?], such that (4.1) holds. One observes very easily that: 

(i) all odd values of M have at least one entry (in bold in Table 1) with g = 

2A( M _ 1 )/2 = (M 2 — l) /4, 5 = 2 M and a = (2 M 2 — l), yielding c = M (M 2 — l) (2 M 2 — l) /2 
with a = M 3 — (3 M — 1) /2; 

(ii) those odd values of a equal to odd integer squares have at least one entry (e.g. 
for M = 17 in Table 1) with g = (a — 1) /8, <5 = (2 a — 1) and a = (2 a + 1), yielding 
c = (a — 1) (4a 2 — l) /8 with M = (y/a — 1) (2a — 1) /2 

(iii) those even values of M equal to twice an integer square have at least one 
entry (in italics in Table 1) with g = M (M 2 — l) /2, 5 = 1, a = y/M/2, yielding 
c = M (M 2 - 1) y/WJ2/2 with a = A M -i. 

These three cases can be generalized respectively to all odd values of M, to all odd 
integer square values of a = (2k + l) 2 , and to all even values of M equal to twice 
integer squares in the following 


Theorem 2. Vfc £ Z + , 3S,a,M,a,c £ Z + such that (2.1) holds : 


(i) if a -- 

= (S 2 - 2) /2, 

with 


(3.1) 

M 

— 

(2k + 1) 

(3.2) 

a 

= 

(2k + l) 3 - (3 k + 1) 

(3.3) 


= 

M* 

2 

(3.4) 

c 

= 

2k (k + 1) (2k + 1) (8k (k - 

(3.5) 


= 

M (M 2 - 1) (2M 2 - 1) 

2 

(ii) if a 

= 5 + 2, with 



(3.6) 


a 

= (2k + l) 2 

(3.7) 


M 

= k (8k (k + 1) + 1) 

(3.8) 



(\/a — 1) (2a — 1) 

2 

(3.9) 


c 

= k(k 2 +ri (m + D‘ 

(3.10) 



(a - 1) (4a 2 - 1) 

8 
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Table 1. First values of M, a, A, S, g, S } a and c for 1 < M < 45 
and 1 < a < 10 5 


a 

A 

s 

9 

23 

72 

578 

2 

25 

135 

735 

15 

96 

490 

9610 

10 

118 

600 

14406 

6 

333 

2352 

112908 

12 

28 

252 

1008 

252 

716 

6480 

518420 

20 

1315 

14520 

1742430 

30 

14 

234 

416 

26 

144 

1950 

22542 

78 

2178 

28392 

4769898 

42 

25 

480 

1080 

120 

3353 

50400 

11289656 

56 

57960 

869505 

3360173145 

105 

9 

289 

361 

1 

120 

2176 

16456 

136 

4888 

83232 

23970888 

72 

153 

2907 

26163 

2907 

680 

12393 

474113 

17 

6831 

129960 

46785690 

90 

14 

504 

686 

14 

144 

3234 

23826 

66 

9230 

194040 

85377710 

110 

12133 

279312 

147476868 

132 

15588 

390000 

243360156 

156 

19643 

530712 

386358518 

182 

81 

2646 

9126 

54 

24346 

706440 

593409810 

210 

29745 

922560 

885657840 

240 

69 

2704 

7396 

4 

133 

4752 

22308 

132 

496 

16368 

261888 

16368 

33 

1617 

2673 

33 

35888 

1184832 

1289097488 

272 

225 

8470 

58870 

70 

42823 

1499400 

1835265906 

306 

50598 

1872792 

2561979798 

342 

111 

5070 

17280 

30 

59261 

2311920 

3514118780 

380 

3276 

131820 

10860720 

780 

68860 

2824080 

4744454820 

420 

64 

3549 

7581 

21 

79443 

3416952 

6314527758 

462 

176 

8910 

39710 

110 

91058 

4098600 

8295566906 

506 


G 

c 

17 

204 

7 

315 

31 

2170 

49 

2940 

97 

16296 

2 

504 

161 

57960 

241 

159060 

4 

312 

17 

6630 

337 

368004 

3 

720 

449 

754320 

5657 

54052635 

19 

323 

11 

5984 

577 

1412496 

3 

8721 

167 

76653 

721 

2465820 

7 

588 

19 

8778 

881 

4070220 

1057 

6418104 

1249 

9742200 

1457 

14319396 

13 

4914 

1681 

20474580 

1921 

28584480 

43 

4472 

13 

10296 

4 

65472 

9 

2079 

2177 

39081504 

29 

22330 

2449 

52457580 

2737 

69267996 

24 

9360 

3041 

90135240 

118 

1196520 

3361 

115752840 

19 

5187 

3697 

146889204 

19 

18810 

4049 

184391460 
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(in) if 5 = 1 and a = k with k > 1 and 


(3.11) 

(3.12) 

(3.13) 

(3.14) 

(3.15) 


M 

a 


c 


2 k 2 

k 2 ( 2k 2 - l) 
M (M — 1) 

2 

k 3 (4A: 4 - 1) 
M (M 2 - 1) 



Proof. For k, 5, a, g , a, M > 1, c £ Z + , g £ Q + , for (2.1) to hold: 

(i) let a = (<5 2 — 2) /2. Replacing in (2.13) yields 

(3.16) g = ^ ^5 4 - 4<5 2 + 4 + ^5 S - 85 6 + 86 4 + 1 6S 2 " 2 ) + 16^ 

For the polynomial in 5 under the square root sign to be a square, let (5/M) 2 — 2 = 
2, i.e. <5 = 2 M, yielding 5® — 85 6 + 85 4 + 32<5 2 + 16 = (5 4 — 4 5 2 — 4) 2 , giving 
g = M 2 (<5 2 — 4) /4<5 2 = (M 2 — l) /4 where 5 was replaced by 2M. As g £ Z + , M 
cannot be even and must be odd, i.e. M = 2k + 1 Vfc £ Z + , yielding g = k ( k + 1), 
5 = 2 M = 2 (2k + 1), a = 2M 2 — 1 = 8fc (k + 1) + 1. Replacing in (2.2) and (2.3) 
yields directly (3.2) to (3.5). 

(ii) Let a = 5 + 2. Replacing in (2.13) yields 

(3.17) 9= M-(m (6 + 2) 2 + \jM 2 ((<5 + 2) 4 - + S^j 

For the polynomial in 5 under the square root sign to be a square, let M = q5 

and 5 = 8 q(q + 1) + 1 with q £ Q+. It yields then ^ jq 2 ^(5 + 2) 4 — + 6 2 = 

5 2 — 4g (5 + 1), giving g = q (q + 1) /2. As g £ Z + , q £ Z + and let q = k, yielding 
g = A*,, 5 = 16Afc + 1, a = 16A& + 3, M = k5 = k (16A*, + 1) =(3.7). Replacing 
in (2.11) and (2.10) yields (3.6) and (3.9). Replacing k in function of a from (3.6) 
yields also (3.8) and (3.10). 

(iii) Let <5 = 1 and a = k with k > 1. Then (2.13) reads 

M (k 2 M + JM 2 (fc 4 - 1) + l) 

(3-18) g = -A--- L 

which takes integer values if M = 2k 2 , yielding g = k 2 (4fc 4 — l) = M (M 2 — l) /2. 
Replacing in (2.11) and (2.10) yields directly (3.12) to (3.15). □ 


The case (i) of Theorem 2 confirms the statement of Stroeker ([37], p. 297) about 
all odd values of M (in bold in Table 1) having a solution to (2.1). The first 50 000 
values of (M, a, c) for this case (i) are given in [28, 29] 

For the case (ii) of Theorem 2, Table 2 gives the first five values of M, a, g, 6, a and 
c. The first 50 000 values of (M, a, c) for this case (ii) are given in [30, 31]. 

As there exist other triplets of values of (M, a, c ) such that (2.1) holds with M = kd 
for the same value of g = A*,, the case (ii) of Theorem 2 can be generalized to other 
values of a as follows. Table 3 shows some of these values for 1 < k < 2 and 
1 < M n , a n < 10 5 and indexed by n for increasing values of M n . 
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Table 2. Values of M,a,g,6,a,c for 1 < k < 5 for case (ii) of 
Theorem 2 


k 

M 

a 

9 

(5 

a 

c 


17 

9 

1 

17 

19 

323 

2 

98 

25 

3 

49 

51 

7497 

3 

291 

49 

6 

97 

99 

57618 

4 

644 

81 

10 

161 

163 

262430 

5 

1205 

121 

15 

241 

241 

878445 


Table 3. Values of M n ,a n ,g,5 n ,a n ,c n for 1 < k < 2 and 1 < 
M n ,a n < 10 5 


Q 

m 

M n 


B 

Sn 


Cn 

D 

D 

17 

9 

i 

17 

19 

323 

D 

O 

305 

153 

i 

305 

341 

104005 

D 

H 

5473 

2737 

i 

5473 

6119 

33489287 

D 

El 

98209 

49105 

i 

98209 

109801 

10783446409 

n 

B 

98 

25 

11 

49 

51 

7497 

m 

B 

4898 

1225 

m 

2449 

2549 

18727503 


It is easily seen from Table 3 that, for each value of k, 

(3.19) <j n S n = <Jn—i V S n —i 

(3.20) a n + 5 n = ^2 (2k + 1) — 1^ cr n —i + ^2 (2k +1) + 1^ S n - i 
with do = d'o = 1, yielding the recurrence relations 

(3.21) 5 n = 2 (2k + l) 2 S n -i — S n -2 

(3.22) a n = 2 (2k + l) 2 a n -i - cr n -2 

i.e. S n and a n fulfill the Diophantine equation (2k (k + 1) + 1) S‘ 2 n — 2k (k + 1) a 
1. The general solutions of this Diophantine equation can be expressed in function 
of Chebyshev polynomials of the second kind U n ^(2 k + l) 2 ^j as 

(3.23) 5 n = U n ((2 k + l) 2 ) - U n —\ [(2k + l) 2 ) 

(3.24) a n = U n [(2k + l) 2 ) + U n - 1 [(2k + l) 2 ) 

These results are generalized for all k and n £ Z + in the next 


3 to 
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Theorem 3. \/k,n € Z + , 35 n , a n , M n , a n , c n £ Z + , M n > 1, such that (2.1) holds 
If &n = + (6% — 1) /2k ( k + 1), with 

(3.25) 

(3.26) 


(3.27) 

(3.28) 

(3.29) 

(3.30) 


M n = kS n 


= k 


u n ((2 k + l) 2 ) - U n _1 ((2 k + l) 2 ) 


CL n — 


$n + 1 


Un ((2 k + l) 2 ) - U n —\ ((2 k + l) 2 ) + 1 


c n — d Tl 


k (k + 1) ((2k (k + 1) + 1) 81 - 1) 


k (k T 1) U2n ^(2/c + 1) ^ 


Proof. For fc,n, S n , a n , g, M n , a n , c n £ Z+, M n > 1, let 
a n = \/5n + (S^ — 1) /2k (k + 1). Replacing in (2.13) yields 


(3.31) 





M2 


/ (2fc (fc + 1) + 1) J2 — 1 

\ 2fc (k + 1) 



For the polynomial in S n under the square root sign to be a square, let M n = kS n , 
yielding after simplification g = k (k + 1)/2 = A*,. Replacing further in (2.2) 
and (2.10) yields (3.27) and (3.29). Replacing further S n by (3.23) yields (3.26) 
and (3.28). Replacing 5 n and a n in (2.10) yields (3.30), noting that S n a n = 

^(2 k + l) 2 ^j — U%_ 1 {(2k + l) 2 ^ = U 2 n ^(2 k + l) 2 ^ as can be shown by replac¬ 
ing U% (x) and Uf _ 1 (x) in function of Chebyshev polynomials of the first kind, 
respectively T 2 n +2 (x) and T 2 n ( x ) (see e.g. [36]) and simplifying appropriately. □ 


For the case (iii) of Theorem 2, the first 50 000 values of (M, a, c) are given in [32], 
There exist also other triplets of values of (M, a, c ) such that (2.1) holds with <5=1 
and a = k. Table 4 shows some of these values for 1 < k < 5 and 1 < M n , a n < 10 5 
and indexed by n for increasing values of M n . The generalization of the case (iii) 
of Theorem 2 to other values of M is included in a following more general theorem. 
This theorem will use the solutions of simple and generalized Pell equations, that 
are recalled in the next section. 


4. Pell equations: A Reminder 
Pell equations of the general form 


(4.1) X 2 - DY 2 = N 

with X, Y, N £ Z and square free D £ Z + , i.e. VD Z + , have been investigated in 
various forms since long (see historical accounts in [7, 14, 38, 13]) and are treated 
in several classical text books (see e.g. [20, 21, 39, 11] and references therein). A 
simple reminder is given here and further details can be found in the references. 
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Table 4. Values of M n ,a n , g n ,c n with 5 = 1 and a = k 
for 1 < k < 5 and 1 < M n . a n < 10 5 


k 

n 

M n 

a n 

<?n 

Cn 

2 

~r 

8 

28 

252 

504 

2 

2 

63 

217 

15624 

31248 

2 

3 

496 

1705 

968440 

1936880 

2 

4 

3905 

13420 

60027660 

120055320 

3 

~r 

18 

153 

2907 

8721 

3 

2 

323 

2737 

936054 

2808162 

3 

3 

5796 

49105 

301406490 

904219470 

4 

~r 

32 

496 

16368 

65472 

4 

2 

1023 

15841 

16728096 

66912384 

5 

~r 

50 

1225 

62475 

312375 

5 

2 

2499 

61201 

156062550 

780312750 


For N = 1, the simple Pell equation reads classically 
(4.2) V 2 - DY 2 = 1 


which has, beside the trivial solution (Xo,Yo) = (1,0), a whole infinite branch of 
solutions Mn £ Z + given by 


(4.3) 

(4.4) 


Xr, 


Y n 


+ v^DVi)” + (V-i - v^DVi) 
2 

(X x + v^DVi)” - (V-i - v^DVi) 
2 s]D 


where (Xi , Y \) is the fundamental solution to (4.2), i.e. the smallest integer solution 
different from the trivial solution (X\ > 1, Y) > 0, £ Z + ). Among the five methods 
listed by Robertson [35] to find the fundamental solution (Xi,Yi), the classical 
method based on the continued fraction expansion of the quadratic irrational y/~D 
introduced by Lagrange [12] is at the core of several other methods. It can be 
summarized as follows. One computes the r th convergent (p r /q r ) of the continued 
fraction [cto;a i,..., ay, a r+ i,...] of y/~D which becomes periodic after the following 
term a. r +i = 2«o if x/~D is a quadratic surd or quadratic irrational (i.e. y/D (j Z + ) 
and with ao = V^dJ i.e. the greatest integer < y/~D . The terms pi and g* can be 
found by the recurrence relations 


(4.5) pi = a.iPi -1 + Pi -2 , Qi = oti.qi-i + qi -2 

with p -2 = 0,p_i = l,<Z -2 = 1,9-1 = 0. The fundamental solution is then 
(Xl,Yi) = (p r ,q r ) if r = 1 (mod 2) or (Xi,Yi) = (p 2r +i, <? 2 r+i) if r = 0 (mod2). 

For the general case of N ^ 1, the generalized Pell equation (4.1) can have either 
no solution at all, or one or several fundamental solutions (Xi, Y\), and all integer 
solutions, if they exist, come on double infinite branches that can be expressed in 
function of the fundamental solution(s) (Xi,Yi) and (Xi, —Yi). Several authors 
(see e.g. [12, 6, 20, 19, 35, 16, 17, 11] and references therein) discussed how to find 
the fundamental solution(s) of the generalized Pell equation, based on Lagrange’s 
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method of continued fractions with various modifications (see e.g. [22]), and further 
how to find additional solutions from the fundamental solution(s). The method 
indicated by Matthews [17] based on an algorithm by Frattini [8, 9, 10] using 
Nagell’s bounds [20, 18], will be used further. 

Noting now ( Xf,yf ) the fundamental solutions of the related simple Pell equation 
(4.3), the other solutions ( X n ,Y n ) can be found from the fundamental solution(s) 
by 

(4.6) X n + VDY n = ± [X r + VDY^ (x f + v^Dy/) " 

for a proper choice of sign ± [35], or by the recurrence relations 


(4.7) 

(4.8) 


X n — XfX n — i + Dy fY n — i 
Y n — XfY n —i H- yfX m -i 


that can also be written as 


(4.9) 

(4.10) 


X n — 2 x f X n _i X n _2 

Y n = ‘lx fY n —\ Y n —2 


or in function of Chebyshev’s polynomials of the first kind T„_i ( Xf ) and of the 
second kind U n -2 ( Xf ) evaluated at Xf (see [23]) 


(4.11) X n = X 1 T n _ 1 {x f ) + DY iyf U n -2 (x f ) 

(4.12) Y n = X 1 y f U n - 2 {x f ) + Y 1 T n - 1 ( Xf ) 

For N = ij 2 an integer square, the generalized Pell equation (4.1) admits always 
integer solutions. The variable change (X',Y') = {{X/if), (Y/rf)) transforms the 
generalized Pell equation in a simple Pell equation X' 2 — DY' 2 = 1 which has integer 
solutions (Xf. Yf). The integer solutions to the generalized Pell equation can then 
be found as (X n ,Y n ) = ( i)X' n ,r)Yf ), or from ((4.9),(4. 10)) with (X 0 ,F 0 ) = (1,0) 
and (Xl,Y i) = ( r)X f ,r)y f ), yielding simply 

(4.13) X n = r]T n (x f ) 

(4.14) Y n = i)yfU n _i (xf) 

(which is also valid for the simple Pell equation (4.2) with ij = 1). Note however that 
not all solutions in (X, Y) may be found in this way (see e.g.[39]) and, depending 
on the value of D, other fundamental solutions may exist. 


5. Step 3: A More General Approach 

The three cases of Theorem 2 can now be generalized as shown in the next theorem, 
that includes also the general method to find all solutions ( M , a, c) such that (2.1) 
holds. 
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Theorem 4. For V<5 ,n £ Z + , 3a,a,M,M n ,c,Xf,yf,Xi,Yi,D,N £ Z + , k £ Q + , 
with gcd (5, a) = 1 ,k = {o'/5) > 1, such as (2.1) holds with 

(5.1) M n = X l y f U n _ 1 {x f )FY l T n {x f ) 

(5.2) a n = -^[(o 2 - 6 2 ) (X 1 +(a 2 + S 2 )Y 1 )y f U n - 1 {x f ) + 

(X! + (o 2 -6 2 )Y 1 )T n (xf)+S 2 ] 

(5.3) c n = ^3 {Xiy f U n -i (x f ) + Y\T n {x f )) 

[(o 2 Xi + (o 4 - 6 4 ) Y ± ) y f U n _i (xf) + (X ± + o 2 Y\) T n (*/)] 

where ( Xf,yf ) and {X\,Y\) are the fundamental solutions of respectively the simple 
(4-2) and generalized Pell equations (4-1) with D = (o 4 — 5 4 ) and N = S 4 , and 
T n ( Xf ) and U n ( Xf ) Chebyshev’s Polynomials of the first and second kind evaluated 
at Xf. 

Proof. Let n,M,M n ,a,a n ,c,c n ,a,5,Xf,yf,Xi,Yi £ Z+, k £ Z*, ti,C,C n £ Q + , 
with M > 1, gcd (5, a) = 1, k = (cr/<5) > 1 . As g (2.13), a (2.2) and c (2.3) in 
Theorem 1 must be integers, the condition for the polynomial( M 2 ( k 4 — 1 ) + 1 ) 
under the square root sign in (2.13), (2.2) and (2.3) to be a squared integer or a 
squared rational allows to find for which values of M (2.1) holds. Let 

(5.4) M 2 (k 4 - 1) + 1 = C 2 
which can be rewritten as a Pell equation as 

(5.5) C 2 - (k 4 - 1) M 2 = 1 
or as 1 

(5.6) (6 2 C) 2 - (o 4 ~6 4 )M 2 =6 4 

which is a generalized Pell equation that always admits at least one fundamental so¬ 
lution as the right hand term is a squared integer. Noting ( Xf,yf ) the fundamental 
solutions of the simple Pell equation (4.2) and {X\, Yj) the fundamental solution(s) 
of the generalized Pell equation (4.1) with X = (S 2 C), Y = M, D = (o 4 — 6 4 ) and 
N = 5 4 , all solutions can be found by (4.6) or ((4.11),(4.12)) Vn £ Z + yielding 

/c ^ (( ct4 ~ <^ 4 ) YiyfU n -i (xf) + XiT n (#/)) 

(O- t) W — ^2 

and (5.1). Then replacing M by M n (5.1) in (2.2) and (2.3) yields directly (5.2) 
and (5.3). □ 

Note that, although (5.1) yields all integer solutions in M to (5.6), some of them 
do not yield integer solutions to a n (5.2) and c n (5.3) and must be rejected. 

For (5 = 1, (5.6) is a simple Pell equation similar to (4.2). It is easy to see that 
its fundamental solution is {xf,yf) — (o 2 , l) (see e.g. [22]). (5.6) admits then an 
infinitude of solutions Vn £ Z + for each integer value of cr, that can be written as 

1 Note that (5.5) could also be written as C 2 — {a 4 — <5 4 ) (M/S 2 ) 2 = 1, but not all solutions 
may be obtained in this way. 
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((4.3), (4.4)) or ((4.13), (4.14)), yielding C n = T n (cr 2 ) and 
M n = U n -i (cr 2 ) 

U n ( O 2 ) - U n -1 (u 2 ) + 1 

8 „ = ----- 

n 2 

crU n (a 2 ) U n -1 (a 2 ) 

Cn ~ 2 

where the relation U n (a 2 ) = T n (cr 2 ) + a 2 U n -i (cr 2 ) (see e.g. [36]) was used in 

(5.9) and (5.10). Note as well that g n = U n (a 2 ) U n -i (cr 2 ) /2 as can be found from 

(2.10) or (2.11). 

This generalizes the case (iii) of Theorem 2 and gives an infinitude of solutions 
(M n , a n , c n ) Vn £ Z + for <5 = 1 and for each value of er, Vcr £ Z + . 

For 5 > 1, three of the fundamental solutions are always (Ai,Yi) = (c> 2 ,0) and 
(<r 2 ,±l), corresponding to respectively ( C,M ) = (1,0) and ((cr 2 /<5 2 ) ,±l). De¬ 
pending on the value of D = (cr 4 — d 4 ), other fundamental solutions may exist. 
All solutions in M can then be found by (5.1) Vn £ Z + for (Ai,Yi) = (<5 2 ,0) or 
(<r 2 ,±l), and Vn £ Z* for other fundamental solutions. Furthermore, solutions 
found for (Xi,Yi) = (cr 2 ,— l) and (<5 2 ,0) yield integer values for M n ,a n and c„, 
while the solutions found for (Xi, Yi) = (cr 2 , l), although yielding integer values of 
M n , do not yield integer values for a n and c„, as can be seen easily from (5.2) and 
(5.3), and must be rejected, although these non-integer values satisfy (2.1). 
Relations (5.7) and (5.1) to (5.3) read 

- for (Xi,Yi) = (cr 2 , -l), C n = k 2 T u ( Xf ) - 6 2 (k 4 - l) VfU n - 1 (xf) and 

<J 2 y f U n -1 (. x f ) - T n (x f ) 

T n ( Xf ) - (cr 2 - <5 2 ) yfU n -i ( xf ) + 1 
2 

crSyfU n -1 (xf) (a 2 y f U n -i (xf) - T n (xf)) 

2 

- for (Xi,Yi) = (S 2 , 0), C n = T n (x f ) and 


(5.11) 

M n 

(5.12) 

a n 

(5.13) 

Cn 


(5.8) 

(5.9) 

(5.10) 


(5.14) 

(5.15) 

(5.16) 


Mn 

a n 

Cn 


S 2 yfU n -i ( x f ) 

T n ( Xf ) + (cr 2 - S 2 ) y f U n _1 (a;/) + 1 
2 

aSyfUn-i ( Xf ) (cr 2 yfUn—i (x f ) + (x/)) 

2 


Note that the solutions (5.11) to (5.13) found for (Xi,Yi) = (cr 2 ,—l) are smaller 
than those (5.14) to (5.16) found for (Xi, Yi) = (<5 2 ,0) for a same value of n > 1. 
Furthermore, if 2 g is a multiple of M, 2 g = gM, y £ Z + , (like for M = 33, 35,42...), 
the fundamental solutions of the simple Pell equation are ( Xf,yf ) = ((/xa 2 — M) , A«), 
as can be easily verified. 
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If g = Am (like for M = 
(5.17) ( x f ,y f ) = 


(5.18) 


5,15,...) or Am-i (like for M = 5,13,15,17,40,...), then 
(((M±1)ct 2 -M),(M±1)) 

O-2-l 


^2 |(5 4 — 2o - 2 + l| 

2 (" 2 -i) v 

|< 5 4 — 2 ct 2 + 1 | J i 


|<5 4 - 2 ct 2 + l| 


- I 


with M = (<5 4 — l) / (1 2 cr 2 — S 4 — l|) (see [33]) with the + (respectively —) sign in 
(5.17) for g = A m (resp. Am-i) and vertical bars denote the absolute value. 
Values of ( Xf,yf ) and (X\,Yi) from [17] yielding solutions (M n , a n , c n ) for the first 
values of M, 1 < M < 45 and 1 < a < 10 5 of Table 1 can be found in [33]. 

For the case (i) of Theorem 2 with M = (2k + 1), g = (M 2 — l) /4, <5 = 2 M and 
a = (2 M 2 - l), one of the other fundamental solutions is always 


(5.19) 


(5.20) 


(*i ,Yi) = 


(5(o- 2 -2)\ J 

2 j ’2 


[M ((2 M 2 - l) 2 - 2^ ,M) 


as can be easily verified. More generally, all solutions for this case can be found by 
(3.1) to (3.5). 

For the case (ii) of Theorem 2 with a = (2k + l) 2 , g = A*,, S = 16A& + 1 = 2a — 1, 
<7 = 16Afc + 3 = 2a + 1, M = kS = k (16A*, + 1 ), one of the other fundamental 
solutions is always known and can be expressed in function of 5, k or a as 


(5.21) (Xr.Vr) = 

(5.22) 

(5.23) 



,2 S 


((4(16A fc + l) ( 8 fc( 2 fc + l ) 3 + l)) ,2(16A fc + l)) 
(( 2 a - 1 ) a 2 (l - a" 1/2 ) , 2 ( 2 a - 1 )) 


6. Recurrent Relations 

The sets of solutions (M n ,a n ,c n ) are obviously not independent. As (5.1) to (5.3) 
are linear combinations of Chebyshev polynomials, one has also the general recur¬ 
rence relations 

(6.1) M n = 2x f M n _i — M n _2 

(6.2) a n = 2xfa n -i - a„_ 2 - (xf - 1) 

(6.3) Cn = 2 ( 2x 2 f - 1 ) C n - 1 - Cn -2 + Sa 3 yj 

among values of M n , a n , c n calculated for the same fundamental solutions (A'i, Yi). 
These relations are immediate from (5.1) to (5.3) and the recurrence and other 
formulas for Chebyshev polynomials (see e.g [36]). For the sake of the recurrence, 
initial values of (M n ,a n ,c n ) for n = 0,1 are shown in Table 5 for the cases (i) 
M„ = kS n of Theorem 3, (ii) <5 = 1 and (xf,yf) = (a 2 , l), (iii) 8 > 1 and (X l5 Yi) = 
(cr 2 ,—l), (iv) 8 > 1 and (X\,Y\) = (<5 2 ,0), and (v) the general case for 8 > 1 and 
with other fundamental solutions (V]\ Y{) of the generalized Pell equation. 
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Table 5. Initial values of (. M n ,a n ,c n ) for n = 0,1 for recurrence relations 


*i,Yi 

n 

(. Mn,an,Cn) for M n = kS n , ( Xf,yf ) = ^(2 k + l) 2 ,0j 

[1,0] 

0 

( 1 , 1 , 1 ) 

1 

Jk (8 k (it + 1) + 1), (2k + if , ( k (k + 1) /2) (4 (2k + l) 4 - l) ) 

Xi,Yi 

n 

(M n ,a„,Cn) for 5 = 1, (xf,yf) = (a 2 , l) 

[1,0] 

0 

(0,1.0) 

1 

(l, <r 2 , cr 3 ) 

*i,Yi 

n 

( M n ,a n ,Cn ) for 5 > 1 

cr 2 , — 1 

0 

(- 1 , 1 , 0 ) 

1 

((a 2 y f - x f ) , (x f - (a 2 - 5 2 ) y f + 1) /2, 5ay f (a 2 y f - x f ) /2) 

<5 2 ,0 

0 

(0,1.0) 

1 

(S 2 y f , ( x f + (cr 2 - 5 2 ) y f + l) /2, 5ay f ( a 2 y f + x f ) /2) 

‘.n* 

0 

(Yf, (X* + (cr 2 - 5 2 ) Yf + 5 2 ) /25 2 , crYj* (X* + cr 2 !^*) /25 3 ) 

1 

((*/ + ( ct2 - 52 ) 2//) x i + (^ - <52 ) ( x f + ( ct2 + 'S 2 ) Vf) Y i + 6 2 ) / 252 , 
cr (y f X* + XfYf) ((x f + a 2 y f ) X* + (a 2 x f + (cr 4 - 5 4 ) y f ) Yf) /25 3 ) 


For <5 = 1, one has also the remarkable recurrence relation 


(6.4) M n = M n — i + 2a ra _i — 1 

among all values of M n Vn £ Z + (see Table 7 further in section 7). One has also a 
similar relation among values of M n and a n for 8 > 1 if only those solutions (5.11), 
(5.12) and (5.14), (5.15) calculated respectively for (Xi,Y\) = (cr 2 ,— l) and (<5 2 ,0) 
are considered. 

If all solutions(M„, a n , c n ) are ordered in increasing value order and indexed by 
a new index j £ Z + , one obtains simply j = n for solutions for <5 = 1, and for 
8 > 1 , if there are no fundamental solutions other than (Xj , Y t ) = (<5 2 ,0) and 
(cr 2 , ±l), one obtains j = 2n — 1 for the solutions (5.11) to (5.13) and j = 2 n for 
the solutions (5.14) to (5.16). If fundamental solutions other than (Xi,Yi) = (<5 2 ,0) 
and (cr 2 ,±l) exist, the solutions ( M n ,a n ,c n ) (5.1) to (5.3) calculated with these 
other fundamental solutions (including for n = 0 ) have to be inserted accordingly 
(see example further in Section 7). 

For 8 > 1, if only those value of ( Mj,aj,Cj ) (5.11) and 5.13) and (Mj + \, aj+i, Cj+i) 
(5.14 to 5.16) calculated respectively with (X\,Yi) = (cr 2 , — l) and (<5 2 ,0) are 
considered, two such sets of solutions are called a “recurrent pair” and obviously, 
for <5 = 1, all sets of ( Mj,aj,Cj ) solutions form “recurrent pairs”. The following 
theorem give other remarkable recurrent relations. 

Theorem 5. For V<5 £ Z + , 3cr, k, j,a,j, Mj,Cj £ Z + with gcd (<5, cr) = 1, k = 
(cf/ 8) > 1, and such as (2.1) holds, if (Mj,aj,Cj) and aj+i, c^+i) fo7~m a 
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“recurrent pair”, then 

(6.5) 

Mj +1 

( 6 . 6 ) 


(6.7) 

Mj 

( 6 . 8 ) 

aj+i + aj 

(6.9) 

a j +1 “ a j 

( 6 . 10 ) 

C J 

( 6 . 11 ) 

Cj +1 + C j 


Mj 2 cij — 1 

K 2 Mj + ^Mf ( K 4 - 1) + 1 = K 2 Mj + Cj 
(2k“ — l) Adj+i — 2Oj+i ~h 1 

k 2 / j _[_ i — Mj -{-1 

(/v 2 - 1) M j+1 

nMjMj + \ 

2 

k 3 M 2 + i 


Proof. Let a,S,j,n,k,a,M,c £ Z + with gcd (<5, cr) = 1, k = (cr/ <5) > 1, and 
( M n ,a n ,c n ) solutions of (2.1), yielding (Mj,aj,Cj) and (Mj + \, ay+i, Cj+i) to be 
a “recurrent pair” with n = j for <5=1, and, for 5 > 1, n = j or n = j + 1 respec¬ 
tively for ( M n ,a n ,c n ) (5.11) to (5.13) for (Xi,Yi) = (<r 2 ,—l) or (5.14) to (5.16) 
for (Xi,Yi) = (<5 2 ,0). Then, 

(i) (6.5) is immediate from (5.8) and (5.9) for <5=1, and from (5.11), (5.12) and 
(5-14) for S > 1. 

(ii) ( 6 . 6 ) is immediate from (2.11), (2.13) and (5.4). 

(iii) For 5 = 1, replacing Mj + 2 by the recurrence relation (6.1) in (6.5) written for 
Mj +2 yields directly (6.7). For S > 1, in (2k 2 — l) M ;+1 — 2a J+ i + 1, replacing 
Mj + 1 and a J+ i by (5.14) and (5.15) yields 


(2k 2 - 1) M j+1 - 2a j+1 + 1 = (2<r 2 — <5 2 ) yfU n -\ (xj) — 

(6.12) (T n (x f ) + (a 2 - S 2 ) y f U n —\ (x f )) 

(6.13) = cr 2 y f U n - 1 ( Xf ) - T n ( x f ) 


which is M n (5.11). Therefore (6.7) holds for 5 > 1. 

(iv) Summing and subtracting 1 and aj extracted respectively from (6.7) and 
(6.5) yield directly ( 6 . 8 ) and (6.9). 

(v) (2.3) and ( 6 . 6 ) yield directly (6.10). 

(vi) Replacing Mj from (6.7) in (6.10) and replacing with Cj +1 from (6.10) yield 

directly ( 6 . 11 ). □ 


This theorem means that once a solution (M„,a„,c n ) has been found for S > 1 
from (5.11) to (5.13) with (Xi,Yi) = (<r 2 ,—l), the other solution ( M n ,a n ,c n ) of 
the "recurrent pair" can be found directly from (6.5) to (6.11) without having to 
be calculated from (5.14) to (5.16) for (Xl, Yi) = (<5 2 ,0). 


7. Summary and Examples 

In summary, to find all solutions (M, a, c ) such that the sum of M consecutive cubed 
integers starting from a > 1 equal a squared integer c 2 , the following approach is 
proposed. 

1) Calculate first all solutions for odd integer values of M = (2k + 1) Vfc £ Z + by 
(3.1) to (3.5). Table 6 shows the first ten values of an infinitude of solutions. 

2) Second, calculate all solutions for the cases M n = kS n given by Theorem 3 
Vfc,n £ Z + either by (3.26), (3.28), (3.30) in function of Chebyshev polynomials, 
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Table 6. Values of (M, a, c ) for M = (2k + 1) , 1 < k < 10 


(3, 23, 204), (5, 118, 2940), (7, 333, 16296), (9, 716, 57960), (11, 1315, 159060), 
(13, 2178, 368004), (15, 3353, 754320), (17, 4888, 1412496), (19, 6831, 2465820), 
(21, 9230, 4070220) 


Table 7. Values of (M n , a n , c n ) for 5 = l,2<a<5 and 2 < n < 6 


cr = 2 

(7 = 3 

(8, 28, 504) 

(63, 217, 31248) 

(496, 1705, 1936880) 

(3905, 13420, 120055320) 

(30744, 105652, 7441492968) 

(18, 153, 8721) 

(323, 2737, 2808162) 

(5796, 49105, 904219470) 

(104005, 881145, 291155861205) 

(1866294, 15811497, 93751283088567) 

a = 4 

(7 = 5 

(32, 496, 65472) 

(1023,15841, 66912384) 

(32704, 506401, 68384391040) 

(1045505, 16188976, 69888780730560) 
(33423456, 517540816, 71426265522241344) 

(50, 1225, 312375) 

(2499, 61201, 780312750) 

(124900, 3058801,1949220937250) 

(6242501, 152878825, 4869153120937875) 
(312000150, 7640882425, 12163142546881874625) 


or by the recurrence relations (6.1) to (6.3) with (. Xf, y /) = ^(2 k + l ) 2 , Oj and the 
initial values of Table 5 (see Table 3). 

3) Third, calculate all solutions for the case S = 1 Vcr £ Z + either by (5.8) to (5.10), 
or by the recurrence relations (6.1) to (6.3) with (a ’f,yf) = (cr 2 , l) and the initial 
values of Table 5. Table 7 shows the first five values of the infinitude of solutions 
for 2 < a < 5. 

4) Fourth, for all other cases with S > 1, find the fundamental solutions of the 
simple and generalized Pell equations with D = (cr 4 — <5 4 ) and N = S 4 for all 
values of cr £ Z + and of <5 £ Z + such that 1 < <5 < cr and gcd (S, a) = 1. 

4.1) If there are no other fundamental solutions than (X\,Y\) = (5 2 , 0) and (cr 2 , ± 1 ). 
then g = fj,M/2 (//, € Z + ) and (. Xf,yf ) = ((/i a 2 — M) , /x); if g = A M or A M _ 1; 
then Mi = (i5 4 — l) / (|<5 4 — 2cr 2 + l|) and (5.18) gives (xf,yj). Then (5.11) to 

(5.16) yield an infinitude of integer solutions (. M n ,a n ,c n ) Vn € Z + with (Vi, Yi) = 
(cr 2 , —l) , (<5 2 ,0). Alternatively, the recurrence relations (6.1) to (6.3) with the ini¬ 
tial values of Table 5 or the recurrence relations of Theorem 5 can be used. Solutions 
have to be ordered then in increasing order of Mj. 

For a = 3 and S = 2, g = 120 = A 15 , D = 65 and N = 16 yielding (xf,yf) = 

(129.16) and only three fundamental solutions (Xi,Y\) = (4,0) and (9, ±1) (see 
[17]). The first ten solutions are shown in Table 8 for (Xi, Yi) = (9, —1) and (4, 0), 
arranged by Mj increasing values and with respectively j = 2n — 1 and j = 2 n. 

4.2) If there are fundamental solutions (X\. Yi) = (X*,±Y*) other than (A 2 ,0) 
and (cr 2 , Al), then ( Xf,yf ) has to be calculated separately (see [17]). All integer 
solutions ( M n ,a n ,c n ) are found by (5.1) to (5.3) for (Xi,Yi) = (A]*',Y 1 *) Vn £ Z* 
and by (5.11) to (5.16) for (Vi, Yi) = (cr 2 ,—l) , (<5 2 ,0) Vn G Z + . Alternatively, the 
recurrence relations (6.1) to (6.3) can be used with the initial values of Table 5. 
Solutions have to be ordered then in increasing order of Mj. 
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Table 8. Values of (. Mj,aj,Cj ) for S = 2, a = 3 with ( Xf,yf ) = 
(129,16) and 1 < n < 5 with j = 2n — 1 and j = 2 n respectively 
for (Vi, Yi) = (9,-1) and (4,0) 


-Yi,Yi 

n 

3 

{Mj , aj , Cj ) 

9, -1 

i 

1 

(15, 25, 720) 

4, 0 

i 

2 

(64, 105, 13104) 

9, -1 

2 

3 

(3871, 6321, 47938464) 

4, 0 

2 

4 

(16512, 26961, 872242272) 

9, -1 

3 

5 

(998703, 1630665, 3190880053872) 

4, 0 

3 

6 

(4260032, 6955705, 58058190109584) 

9, -1 

4 

7 

(257661503, 420705121, 212391358097903424) 

4, 0 

4 

8 

(1099071744, 1794544801, 3864469249201901760) 

9, -1 

5 

9 

(66475669071, 108540290425, 14137193574521767668240) 

4, 0 

5 

10 

(283556249920, 462985602825, 257226802107318794853360) 


Table 9. Values of ( Mj,a,j,Cj ) for 5 = 3, a = 4 with ( Xf,yf ) = 
(2024,153) and 0 < n < 3 with j = 3n + 1, j = 3n — 1 and j = 3 n 
respectively for (Xi,Yi) = (159,12), (16, —1) and (9,0) 


-YiTi 

n 

3 

(Mj,aj,Cj) 

159,12 

0 

1 

(12, 14, 312) 

16,-1 

1 

2 

(424, 477, 389232) 

9,0 

1 

3 

(1377, 1548, 4105296) 

159,12 

1 

4 

(48615, 54635, 5117020440) 

16,-1 

2 

5 

(1716353, 1928872, 6378077594592) 

9,0 

2 

6 

(5574096, 6264280, 67270624549920) 

159,12 

2 

7 

(196793508, 221160443, 83849042274507096) 

16,-1 

3 

8 

(6947796520, 7808071356, 104513105644422184080) 

9,0 

3 

9 

(22563939231, 25357801869, 1102316769603603585072) 

159,12 

3 

10 

(796620071769, 895257416606, 1373975729120835063673080) 


For a = 4 and S = 3, D = 175 and N = 81 yielding ( Xf,yf) = (2024,153) and 
five fundamental solutions (X\, Yi) = (9,0), (16, ±1) and (159. ±12) (see [17]). The 
first ten solutions are shown in Table 9 for (Vi, Yi) = (159,12), (16, —1) and (9,0), 
arranged by Mj increasing values and with respectively j = 3n ± 1, j = 3?r — 1 and 
j = 3 n. 


8. Conclusions 

The approach proposed in this paper to find all solutions ( M,a,c ) to (1.2) by 
investigating two other parameters, S and a, instead of investigating each individual 
values of M one by one in an increasing order of M values allows to find quite simple 
and elegant general solutions (M, a, c) based on solutions of simple and generalized 
Pell equations involving Chebyshev polynomials. Alternatively, recurrence relations 
can be used in order to simplify the computational part. This approach allows to 
find all possible solutions (M, a, c) to (1.2) but with the drawback that solutions 
are not ordered by increasing M values. 
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However, it is found that there are always at least one solution for every cases of 
all odd values of M, of all odd integer square values of a, and of all even values of 
M equal to twice an integer square. 
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